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We evaluate systematically some new contributions of the QCD scalar mesons, including radiative decay- 
productions, not considered with a better attention until now in the evaluation of the hadronic contributions 
to the muon anomaly. The sum of the scalar contributions to be added to the existing Standard Model predic- 
tions af M are estimated in units 10 -10 to be = 1.0(0.6) [TH based] and 13(11) [PDG based] , where the errors 
are dominated by the ones from the experimental widths of these scalar mesons. PDG based results suggest that 
the value of a^ M and its errors might have been underestimated in previous works. The inclusion of these new 
effects leads to a perfect agreement (< l.lcr) of the measured value a™ p and a^ M from r-decay and implies a 
(1.5 ~ 3.3) a discrepancy between a e ^ v and a^ M from e + e~ — > hadrons data. More refined unbiased estimates of 
a^ M require improved measurements of the scalar meson masses and widths. The impact of our results to a^ M 
is summarized in the conclusions. 



1. INTRODUCTION 

Improved accurate experimental measurement pQ 
and recent theoretical estimates of the muon 
anomaly |2I3I4I5| are now available. The theoreti- 
cal accuracy in |2I3| 1 is mainly attributed to the 
use of the new CMD-2 (Sj around the p mass and 
on BES data [7] around the J/ty region. The im- 
pact of the former data on the determination of 
the muon anomaly is intuitively more important due 
to the low-energy dominance of the anomaly ker- 
nel function |8I9| . However, what is more intrigu- 
ing within this increasing precision is the discrep- 
ancy between the results from e + e~ and r-decay 
data |2l,'{j . which was not the case in the previous 
determinations using preliminary data |4I5| . In this 
short note, we study some other sources of contri- 
butions, not considered with a better attention until 
now, from the scalar mesons. These scalar (qq, glu- 
onia,...) are conceptually fundamental consequences 
of QCD. Though their existence is not precisely con- 
firmed, there are increasing evidences of their find- 
ings in different e + e~ and hadronic experiments 10 . 



1 Though agreeing in the total sum, the results from these two 
estimates differ in each energy region. 



2. ISOSCALAR SCALAR MESONS 

The interest in these I — scalar mesons are 
that they cannot be obtained from usual ChPT 
approaches. They are related to the QCD scale 
anomaly: 

0» = \f3(a s )G 2 + £[1 + 7 m(a s )WiV>i , (1) 

i 

where G" is the gluon field strengths, ipi is the 
quark field; (3(a s ) and "f m (a s ) are respectively the 
QCD /3-function and quark mass-anomalous dimen- 
sion. In this case, arguments based on SU{2) sym- 
metry or its violation used to estimate some pro- 
cesses like e.g. radiative processes cannot be applied 
[HI- Using QCD spectral sum rules (QSSR) fI21T3) 
and low-energy theorems (LET) for estimating the 
mass and its width, it has been shown in |14I15) that 
the wide er(0.60) meson is the best candidate me- 
son (gluonium) associated to the previous interpo- 
lating current (see also JHl)- Its mass is due to the 
gluon component and its large width into 7T7t is due 
to a large violation of the OZI rule 2 (analogue of 
the r/-meson of the U(1)a sector J7|). The oberved 

2 This feature also implies that lattice calculations of the gluo- 
nium mass in pure Yang-Mills or quenched approximation are 
bad approximations and might miss this "unusual" glueball. 
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<r(0.60) and /o(0.98) can be explained by a "max- 
imal quarkonium-gluonium mixing scheme" |14I15| . 
In the following, we study the a and other scalar 
mesons contributions to the muon anomaly. 



3. THE e+e 



S 7 PROCESSES 



Using analytic properties of the photon propagator, 
the general contribution to the muon anomaly from 
the process: 



e e — ► 7* — > hadrons 

can be written in a closed form as 

1 r°° 

4 ad (l.o) = — dt Ki{t) a H (t) , 
4tt j J Am 2 

where: 

• Ki(t> 0) is the QED kernel function 



(2) 
(3) 



dx- 



x 2 {l 



, x- + (t/mf) (1-x)' 
with the analytic form: 

v2 N 



(4) 



Ki(t>Amf) 



z 



1 



z 



f )+(! + «)' x 



log(l + zj) - zi 

+ zi 



1 



zi 



z l log zi, (5) 



with: 



VI = 



t 



4m 2 ' 



zi 



(1 " vi) 



and vi 



1 



4m? 
t 



(6) 



is a monotonically decreasing function of t. For 
large t, it behaves as: 

™2 



i^(i > mf) ~ , 



(7) 



which will be an useful approximation for the anal- 
ysis in the large t regime. Such properties then em- 
phasize the importance of the low-energy contribu- 
tion to af ad (l.o) (I = e, /i), where the QCD analytic 
calculations cannot be applied. 

• c.ff(£) = <r(e + e~ — * hadrons) is the e + e~ — > 
hadrons total cross-section which can be related to 
the hadronic two-point spectral function lmH(t) em 
through the optical theorem: 

er(e + e~ — > hadrons) 



Re+e- 

where 
a(e + e 
Here. 



cr(e + e — > ) 



47ra 2 



MM j = 



3£ 



127rlmn(i) em , (8) 



(9) 



i /d^ (OITJ^^) GC^))* 10) 

<f<f)iW<7 2 ) (io) 



(<T<? 2 



is the correlator built from the local electromagnetic 
current: 



■4m (z) = \u^u - ±drfd ~ \*f8 + 



(11) 



Using Vector Meson Dominance Model (VDM) in a 
Narrow Width Approximation (NWA), one obtains 



,VDM 



(1.0) 



-K 



Ml 



r 



M- 



(12) 



The V — + ■yX (X = a,fo) coupling has been esti- 
mated from the V 77 one given in |14I15| and in the 
PDG data [H| and using the VDM relation: 



gx-yv 



\/2 7 v 



-9x 



77 ' 



(13) 



where e is the electric charge and we use the normal- 
ization 7 P ~ 2.51 ± 0.02. For the scalar mesons, we 
use: 

M a ~ (0.6 - 0.8) GeV , 
T a ^ ~ (0.2 - 0.3) keV QSSR [nUTS] 
~ (3.8 ± 1.5) keV PDG QHH] , 

r /o(1 . 4HrT ~ (0.7~5.)keV, (14) 

The ones of the other w and <p radiative widths come 
from PDG [H]. We also use in MeV units [15]: 

i 0.08 x 10~ 3 , T tot ~ 174(59) 
; 0.44 x 10~ 3 , T tot ^ 310(60) 
; 0.48 x 10~ 3 , r tot ~ 150(50) , 



r 



w(1.42)->e+e- 
■ p(1.45)->e+e" 



F 



0(1.68)- 



(15) 



We deduce the results in Tabled where the first set 
of values corresponds to r CT ^ 77 from QSSR and the 
second set of values to the a width from PDG. It is 
clear from Tabled that the results are very sensitive 
to the value of the a mass and 77 width. Though, in 
the "maximal quark-gluonium mixing scheme" of the 
a meson, one favours a small 77 width |14I15| . the 
PDG data |19I20| and some other QCD models still 
allow higher values. For a conservative and unbiased 
estimate, we translate the total sum in Tabled into: 



xlO 10 = 



0.45(0.13) QSSR 
7.30(5.60) PDG . 



(16) 



In order to include this result into the one in and 
for avoiding an eventual double counting, we "de- 
rescale" their result on the top of the u) by 0.888 and 
on the 4> by the factor 0.984 used there for taking into 
account the missing modes. In this way, we obtain 
the corresponding "pure" hadronic decay: 

a^uj + <j) -> 3tt + KK) = 67.04(1.50) x 10~ 10 , (17) 

Using the to and <fi into 777 and ir°j radiative decay 
widths from |19| . we deduce: 

a^uj + <j) -> 777 + tt°7) = 4.36(0.16) x 10~ 10 . (18) 
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Table 1 



scalar+7 contributions to a^ ad *K 



Processes 



a had x 1Q 10 







QSSR 




PDG 






- 0. 




- 0.05 


to —>■ 0-7 




- 0. 




- 0.07 


4> — ► o"7 




0.03 ~ 


0.15 


1. - 10. 


/ (.98) 7 


PDG 




0.01 ±0.00 


^ — > ao(.98)7 


PDG 




0.03 ±0.01 


p(1.45) - 


* (77 


0.01 ~ 


0.02 


0.06-0.42 


w(1.42) - 


-> (77 


0.10 ~ 


0.17 


0.1 - 0.7 


0(1.68) - 


■» (77 


0.14- 


0.20 


0.17- 1.17 


0(1.68) - 


- /o (1-4)7 








Total 




0.45 ±0.13 


7.3 ±5.6 



NEW CONTRIBUTIONS 
SCALAR MESON TO a„ 



OF THE 



Here, we study the effect of scalar mesons via a 
Higgs-like triangle diagram shown in Fig. ^ This 
new hadronic contribution to a^, ad cannot be esti- 
mated using the usual electromagnetic spectral func- 
tion associated to one virtual photon. The evaluation 




*' For completeness, we also include the contri- 
bution of the isovector ao(980). 



Figure 1. Scalar meson contribution to a;. 



Taking into account these different effects, we con- 
clude from this analysis that the inclusion of these 
new radiative decays increases the results in [2] by: 



Aa*~* x 10 



10 



= 0.10(0.13) 
= 6.95(5.60) 



QSSR 
PDG . 



(19) 



PDG data for the widths of the a meson may indicate 
that the scaling factors used in |2| are not sufficient 
for taking into account the different missing modes of 
the lo and <f> mesons. Adding our result in Eq. i|19[l . 
to the recent estimate of a^ ad (l.o) in |2I3I4I5| . the 
lowest order hadronic contributions to becomes 3 
in units of 10" 10 : 



Aiad 



(l.o) = 706.4(6.8) (0. 1) [r- QSSR] , 

= 713.3(6.8)(5.6) [r - PDG] , 

= 683.7(6.4)(0.1) [e + e"-QSSR] , 

= 690.6(6.4)(5.6) [e + e~ — PDG] , (20) 



where the first error is the ones from |2l3l4] . and the 
second one comes from the present analysis. Our re- 
sults, especially the one based on PDG data, suggest 
that the estimate of a^ ad (Lo) and its errors might 
have been underestimated in the previous determi- 
nations. 



3 For the r-decay data, we use the average of the results (in 
units of lCT 10 ) 709.0(5.9) H and 703.6(7.6) E|. For the 
e + e _ data, we use the average of the most recent estimates 
684.0(6.5) |21 and 683.1(6.2) which are however lower by 
about 1.4cr than the previous estimates in 14151 based on pre- 
liminary data. 



of this diagram gives the contribution (see e.g. |21p: 



\9F\' 

16tt 2 



d.r 



x 2 {2-x) 



+ (Ml/ml) (1 - x) 



(21) 



where Ms is the scalar mass and gF is the "effec- 
tive" Sl + l~~ coupling. However, the coupling of these 
scalar mesons to e + e~ is not known, and in the case 
of the isoscalar, the naive argument based on chiral 
symmetry (Higgs-type coupling) may not be applied 
due to the presence of gluon fields in the U(l)v dila- 
ton current given in Eq. QJ. Note that these scalar 
mesons might also be produced in the s-channel in 
e + e~ experiments, most probably via two-photon ex- 
change, which should differ from the usual spectral 
function of a one photon exchange used to estimate 
a^ ad (Lo), and therefore prevent from some eventual 
double counting of the e + e~ — * hadrons data. How- 
ever, we are not also aware of existing analyses of 
the angular distribution which can differentiate a 
scalar from a vector meson in the low-energy e + e~ 
region concerned here (22], while present LEP lim- 
its on electroweak scalar particles (sleptons, squarks) 
given in [211] may not be applied here.. If one esti- 
mates this coupling from the experimental bound on 
the S — > e + e~ width 4 , and use the positiv- 
ity of the contribution, one gets the result in Table 
121 where an universal coupling to l + l~ (I = e,n) 
has been assumed 5 . In order to see the strength of 
this experimental bound, one may assume that the 



4 Note that the difference a% xp - of M = (34.9 ± 28.4) X 10" 12 
from the electron anomaly 1241 leads to a weaker upper bound. 
5 The bound is weaker for a Higgs-type coupling. 
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Table 2 

Contribution of Fig. to a M using the upper bound 
on the leptonic widths compiled by PDG |19| . 





Bp 


x 10 10 


Ts^e+e- [eV] 


<r(.60) 


0. ■ 


- 8.9 


< 20 *) 


/o(-98) 


0. ■ 


- 1.1 


< 8.4 


a (.98) 


0. ■ 


- 0.2 


< 1-5 


/o(L37) 


- 


- 1.1 


< 20 


Total 


0. • 


- 11.4 





*J Due to the non-available experimental lep- 
tonic width of the cr and the unclear separation 
between the a and the /o(1.37), one can safely 
assume that the bounds for the cr and for the 
/o(1.37) are about the same. 



coupling of the scalar to is dominated by the 
one through two photons. Therefore, an alternative 
rough estimate of this effect can be obtained by re- 
lating it to the light-by-light (LL) scattering diagram 
contribution 1251. where we obtain 6 : 



{LL) 



(0.1 - 1.0) x 10" 



(22) 



where we have considered the uncertainties to be 
about one order of magnitude in order to have a 
conservative estimate. This result indicates that the 
present experimental upper bound in Table |2] might 
be very weak, and needs to be improved. We trans- 
late the results given in Eqs. (|22|l and Table [21 into 
(in units of 10~10): 

n Sll 



0.55(0.45) TH 
5.7(5.7) PDG. 



(23) 



Adding the results in Eqs. (|19|l and (|23fl . one finally 
deduces the sum of the additionnal contributions due 
to scalar mesons in units of 10~ 10 : 



1.0(0.6) TH 
13(11) PDG 



(24) 



to be added to the SM predictions a„ . One should 
notice that the present theoretical estimate based on 
QCD spectral sum rules (QSSR) differs by about an 
order magnitude to the one based on the PDG data 
for the (T77 an d scalar meson leptonic widths. The 
QSSR estimate of the 0-77 width is based on the 
picture where the observed a meson emerges from a 
maximal mixing between a gluonium and qq states, 
which then implies a much smaller 77 width. On the 
other, the theoretical estimate of the leptonic width 
is based on an intermediate 77 coupling of the scalar 

6 We use the total LL contribution in order to take into account 
other mesons contributions 



meson to the lepton pair explaining again the relative 
suppression of a such width. A progress in improving 
the accuracy of the scalar meson contributions needs 
solely improved measurements of the scalar mesons 
77 and leptonic widths. In addition, such a program 
is also necessary to clarify the exact nature of the 
scalar mesons which, at present, is, experimentally, 
poorly known. 

5. CONCLUSIONS 

The inclusion of the scalar meson contribution in 
Eq. (|19|) modifies the recent estimate of a^ ad (Lo) 
from |2I3I4| into the one in Eq. (J2UJ), while the di- 
rect exchange of the scalar meson shown in Fig. ^ 
gives the new contribution in Eq. (|23|l . The total 
sum of the scalar contributions is given in Eq. I|24l) . 
Adding our result in Eqs. i|2"U|) and J23J), we deduce 
in units of 10 -10 : 



Jiad 



had 



{1.0) 



707.4(6.8) (0.6) [r - TH] , 
719.4(6. 8)(11) [r — PDG] , 
684.6(6.4) (0.6) [e + e~ - TH] 
696.6(6.4) (11) [eV-PDG] 



(25) 



where the first error is the ones from |2I3I4| . and the 

second one comes from the present analysis. Our 
results (especially the one using PDG data) suggest 
that the value and errors of the hadronic contribu- 
tions to a^ M might have been underestimated in the 
previous determinations. From the above results, we 
deduce our final estimate of the muon anomaly in 
units of 10- 10 : 

al M = 11 659 191.6(6.8) [r - TH] , 
= 11 659 203.6(12.9) [r - PDG] , 
= 11 659 169.2(6.4) [e + e~ - TH] , 
= 11 659 181.2(12.7) [e+e~ - PDG] 

compared with the recent data £Q: 
a ex P = n 65Q 204.0(8.6) x 10" 10 . 
Then, we deduce in units of 10~ 10 : 

= 12.4(11.0) [r-TH] , 

= 0.4(15.5) [t-PDG] , 

= 34.8(10.7) [e+e- - TH] , 

= 22.8(15.3) [e+e- - PDG] . 



(26) 



(27) 



exp 



,S M 



(28) 



One can see that the inclusion of these new effects 
due to scalar mesons improves the agreement be- 
tween a e * p and a^ M : though a such agreement is 
quite good from r-decay (less than 1.1 cr), there 
still remains some disagrement (1.5 to 3.3 a) be- 



tween a 



SM 



and af, xp . More refined esti- 



mates of a^ M require improved measurements of the 
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masses and widths of the scalar mesons, which are, 
at present, the major obstacles for reaching a high- 
precision value of a^ M . These new measurements 
are also necessary for making progresses in our QCD 
understanding of the nature of scalar mesons, which 
play a vital role on our understanding of the symme- 
try breaking in QCD. In addition to some problems 
mentioned in pj] (absolute normalization of the cross- 
section, radiative corrections, SU(2) breakings,...), 
the discrepancy between the r-decay and e + e~ data 
should stimulate improvements of the data and ex- 
perimental searches for new and presumably tiny ef- 
fects not accounted for until now. We plan to come 
back to this point in a future work. 
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